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1*  In  two  previous  memoranda.,  (references  (a)  and  (b))^ methods  were 
described  for  calculating  the  sound  pressure  on  circular  flexural  disks 
and  obtaining  their  radiation  impedance.  Unfortunately,  it  was  only  feas- 
ible to  obtain  the  sound  pressure  by  means  of  an  electronic  computer,  and 
the  radiation  impedance  data  had  to  be  obtained  by  numerically  integrating 
this  computer  data  for  the  sound  pressure.  In  this  memorandum  the  double 
integrals  needed  to  compute  the  radiation  impedance  are  evaluated.  Simple 
expressions  are  given  for  radiation  impedance  suitable  for  direct  hand  cal- 
culation. 

According  to  reference  (a),  radiation  Impedance  for  circular  flexural 
disks,  referred  to  their  average  normal  surface  velocity  u,  is  given  as 


Most  flexural  disk  problems  of  interest  can  be  formulated  using  only  even 
powers  of  r/a.  Using  the  relation 

P(r)=  pcU  (3) 
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The  pressure  components  are  given  in  reference  (a)  for  the  case  when  k is 
0,  2,  and  4. 


Interchanging  integration  and  summation,  ve  have 

* 2(^/2  S £&(*,.)  *,HJy  (6) 


and 


*/rJ  -- 


die  <*Ji 


K~  A 

From  Eou.  (2)  ve  see  that 

1-1=  XX  S.  >2,  (SkY 


1 


. v S"  I 
X d Y 


(7) 


(8) 


n 


and 


l*|:»  1 I"1  rirU  = zU 


^ cLh 
ntl 


We  also  define  the  complex  quantity 


Qiri  * Sko  t 1 1 Kit  2 51 


so  that 


^^4  = 0A5  _*a_f  ^ 2 **K~  *M  Q 

•r-  t-'" 


kA 


(9) 


It  now  remains  to  evaluate  the  following  two  integrals: 

SK„*  J Tb‘  . (10) 


(11) 


(12) 


When  k and  j?  are  even,  we  shall  be  able  to  integrate  S^j  and  Tj^  . Two 
examples  will  be  carried  out,  Qpo  and  Q20  Tables  I and  II  show  and 

Tj^  for  k and  S » 0,  2 and  4.  It  will  be  shown  in  the  Appendix  that 

Qk*  ~ Q n- 

Referring  to  Eou.  (ll). 
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^ 1 y-  **  = L f'"  ±Le',kRoJ*)*<i*  (13) 

“ Z io  £,  e~,kr^<J*  . 


In  reference  (a),  RQ  was  defined  as 
R.o  - rCt>sr  ot  4-  (o-L  - 


<4 


(i*0 


so  that 

- -y.  a os  jc  + (l  — 


Vi 


a 


Therefore. 


(15) 


. , (*  s'tt  -;kk  (•**■*«- + *1 -+.**•»**.}  , 1 

Q"~*  T “nf  J.  Jo  « * 


(16) 


In  the  x-  oc  plane,  the  integration  is  over  a semicircle  with  unit  rad 
ius.  We  now  make  a change  of  variables 

iX  - y <Jir\  oC  j V-  cosol  (17) 

and  the  differential  area  element  x Jx.  dx  becomes  du  dv,  and  ve  have 


When  the  integration  over  v is  performed,  we  have. 


L9) 
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Letting  u = sin  0,  the  last  integral  in  Equ.  (19)  becomes 

/ ' * 


37  (ZhTa.) 


This  last  integration  vas  done  with  the  help  of  one  integration  by  parts,  and 
reference  (c).  S,  (2ka)  denotes  the  Struve  function  of  order  one  and  argu- 
ment 2ka.  Therefore,  we  have 


so  that  finally  we  get 


If  Uo  = 1,  and  all  other  * 0,  then  we  have  simply  the  piston  with  uni 

form  surface  velocity.  Referring  to  Equ.  (12), 

fc/'  a - 7 Cn’  - I—  ^ (2KZx>  /**..» 


These  are  the  familiar  equations  for  the  radiation  resistance  and  reactance 
of  a circular  piston. 


The  same  methods  that  were  used  to  evaluate  Q0o  wiH  be  used  in  obtain 


Referring  to  reference  (a)  for  Zp,  and  again  setting  x = r/a 
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By  again  letting  u = s in  0,  expressing  cos20  as  l-sin20,  and  using  reference 
(c),  we  obtain 


I I ...  t — u 

! 1 — • .i  *• 


J-  - 


*7*.  -i.KcOes 


= 1 r<2te.)  _ S.OCQ)] 

L K Ztav  /J}' 


Finally, 


c » _ 

■J2C-*  X 


_j _ X-UWQ  _ (fc;  f 

Z(K^)  S ~~  (35) 


+ (far-*)  S,(2*Z-) 

(/<«•  J ^ 2-  ('*<'“-)  3 


Considerable  arithmetical  care  must  be  taken  in  using  the  formulae  in 
Table  I.  The  numerical  computation  of  the  and  Tj^  often  involves  the 

subtraction  of  large  comparable  numbers,  which  leads  to  a reduction  in  sig- 
nificant figures  in  the  result.  For  this  reason,  a very  accurate  Table  of 
Bessel  and  Struve  functions  is  required,  with  perhaps  seven  (7)  figures, 
such  as  is  found  in  reference  (d). 

Unfortunately,  this  work  has  not  led  to  a formula  not  involving  inte- 
grals for  the  pressure  distribution  on  a circular  piston  or  on  a flexural 
disk.  However,  it  has  led  to  such  a formula  for  the  integral  of  this  pres- 
sure over  the  disk. 


_____  _ — 7 
J V t ■ 


DAVID  T.  PORTER 
Mathematician 
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APPENDIX 

THE  RECIPROCITY  OF  AND  Q k 

We  wish  to  show  that  Qk||  = Q J?k*  To  do  this,  we  first  write  Q]^  as 


Q**  = \0Z<  * T*  |0  Zfcr  Jr  • 


Integrating  with  respect  to  an  angle  oc  around  a circle, 
/>  , f'  cl  _ f^/ \Q 

bjk*  - 


t 


ITT#-1 


(Al) 


(A2) 


This  double  integral  can  be  considered  as  the  integral  of  Zk(r/a)  over 
the  circular  area  S j/  . Therefore, 


"2tt  ft.1  J Jy 

and  so,  we  wish  to  show  that 

j = y 2^  cr^)K  j s* 

£j?  Sic 


(A3) 


(AA) 


Consider  now  the  mutual  radiation  impedance  between  two  superimposed 
circular  pistons,  k and  $ , which  have  normal  surface  velocities 

V„K  = UK  (A5) 


and 


Viol  - D i e' 


(A6) 


The  mutual  radiation  impedance  between  the  two  pistons,  referred  to  their 
edge  velocities,  Uk  and  Ufl  , will  be,  following  reference  (e),  Equ.  (8>, 


J t,w  Ks 


Because  the  reciprocity  of  Zkj  and  Z ^ k applies  here, 

Zkfl  = Zqk 


(A7) 


(a8) 


i M V'ti  fo)  r*-^~ r f 8^Ktfe)^.(A9) 

vf  Vi  V? 


and 
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Here 


* C$  « 

/*  - t^e'^ 

(A10) 

and 

fy  : Uji  e'^'Zji  J 

Pfcr  - Mhr  4 fc  2Trc- 

(All) 

Therefore,  we  have 

\ jOC 

U*  Vj  ~>sf 

(A12) 

1 

C pc  Oj,  ^ <15tr  . 

‘'SW 

Up  C/fc 

And  so. 


which  was  needed  to  be  shown  to  prove  ^ 


<Mk* 


(A13) 
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TABLE  I 
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TABLE  II 
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